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O ' 

^ , Abstract. 77 Carinae is considered to be a massive colliding wind binary system with 

^ ■ a highly eccentric (e ~ 0.9), 5.54-yr orbit. However, the companion star continues 

' to evade direct detection as the primary dwarfs its emission at most wavelengths. Us- 

ing three-dimensional (3-D) SPH simulations of Car's colliding winds and radiative 
transfer codes, we are able to compute synthetic observables across multiple wavebands 
I for comparison to the observations. The models show that the presence of a companion 

. star has a profound influence on the observed HST/STIS UV spectrum and H-alpha line 

0^ ' profiles, as well as the ground-based photometric monitoring. Here, we focus on the 

\ Bore Hole effect, wherein the fast wind from the hot secondary star carves a cavity in 



(N 



the dense primary wind, allowing increased escape of radiation from the hotter/deeper 
layers of the primary's extended wind photosphere. The results have important impli- 
cations for interpretations of 77 Car's observables at multiple wavelengths. 



■ 1. Introduction 



The 5.54-year periodicity observed in numerous spectral lines, as w ell as the IR, Vi- 
sual, and X-ra y fluxes, s trongly suggests 77 Ca rinae is a binary s ystem (Whiteloc k et al.l 
12004; DamineU et al.ll2008a.b.: ,Fernandez-Lajus et al.ll2010l : ICorcoraa.2011.) . Recently, 
lOkazaki et al. (2008, hereafter O08) modeled 77 Car's RXTE X-ray light curye using 
a three-dimensional (3-D) Smoothed Particle Hydrodynamics (SPH) simulation of the 
binary wind-wind collision. A key point of O08 is that the fast wind of the secondary 
star, r]B, carves a cavity in the dense wind of the primary, 77^, allowing X-rays that 
would otherwise be absorbed to escape into our line-of-sight. This immediately sug- 
gests that if the primary wind is sufficiently optically thick in the UV, Optical, or IR 
waveband, the low-density secondary wind may likewise carve or "bore" a cavity or 
"hole" in the associated wind photosphere, allowing increased escape of radiation from 
the hotter/deeper layers. 
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2. The Bore-Hole Effect 

Any bore-hole effect should depend on (1) how close the wind cavity carved by t]b 
gets to r]A and (2) the apparent size of tja's photosphere. For a bore-hole effect to 
occur, the optically thick primary wind photosphere must extend at least as far as the 
separation distance between the primary and the apex of the wind-wind interaction 
region (WWIR), i.e. if at some point (2) > (1), there is a bore hole. 

The separation between ija and the apex of the WWIR is determined by the dy- 
namics of the colliding winds and the orbit. The absolute minimum distance that t/^'s 
photosphere must extend is the separation between tja and WWIR apex at periastron: 

_ a{l - e) 

^rain — ~, ~~r~i p > (1) 

1 + 1/ 

where e is the eccentricity (0.9), a is the orbital semi-major axis length (15.4 AU), and 
T] is the momentum ratio of tja'^ wind to ijb's wind. For the parameters of O08, rj ~ 4.2 
and j?min ~ 1 AU. Using the primary mass-loss rate M = lO^^M© yr"' suggested by 



Hillier et al.l (l200lL i2006. hereafter HOI and H06), rj ~ 16.67 and /?mm « 1-24 AU. 

The apparent size of /y^'s wind photosphere is determined by radiative transfer 
effects. While there are many ways to define this apparent size, a simple way is to use 
the radial photospheric radius /?phot at which the radial optical depth t = I: 

^Phot = ^ > (2) 

where k is the opacity (units of cm^ g~^, assumed constant for simplicity), M is the 
mass-loss rate of ija, and v is the terminal speed of 77^'s wind. Assuming k = I cm^ g'^ 
and mass-loss rates of M ^ 2.5 x 10""^ Mq yr"' and M ^ IO'^Mq yr"' , /?phot » 1.67 AU 
and 6.71 AU, respectively. Thus, even for low values of k of order unity, a bore-hole 
effect should occur in the t] Cai system (at least at periastron). 

In reality, any bore-hole effect in 77 Car will be wavelength dependent and occur 
over a relatively short period of time around periastron. The wavelength dependence 
is due to the opacity, which determines the size of tja's photosphere. Figure [T] illus- 
trates this, showing an overlaid plot of the normalized intensity as a function of impact 
parameter from the center of tja in each of the wavebands B, V, R, I, J, H, K, L. 

Since the separation between tja and the apex of the WWIR varies with orbital 
phase 0, the bore-hole effect is also orbital-phase dependent. At apastron, the separa- 
tion is ~ 23.6 AU, much larger than tja in either the K- or B-band. Because of the large 
orbital eccentricity, this separation remains large for most of the orbit, > 10 AU for (p 
between 0.5 and 0.925. However, between (p w 0.925 and 1.0 (periastron), the sepai^a- 
tion drops rapidly, from ~ 10 AU to 1.24 AU. At (p ~ 0.97, the separation is ~ 5 AU, 
which is about the same size as tja in the ^-band. Therefore, for the wavebands con- 
sidered here, a bore-hole effect should occur in rj Car at phases 0.97 < < 1.03. 

Because of the different photospheric radii in each waveband, the wind cavity will 
bore different relative distances into tja at a given in a particular waveband. Thus, at 
(p = 0.98, one would expect to see a significant bore-hole effect in the A'-band, but only 
a relatively minor effect, if any, in the B-band. Moreover, because of the differences 
in relative flux probed between wavebands, one might also expect different relative 
changes in brightness before, during, and after periastron. 
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Figure 1. Curves showing normalized intensity vs. impact parameter in each of 
the wavebands (from innermost to outermost) B, V, R, I, J, H, K, L, computed using 
the HOI, H06 CMFGEN model of 774. 




Figure 2. 3-D renderings of the bore-hole effect at orbital phases of apastron (left), 
5 days before periastron (middle), and 155 days after periastron (right). The x and 
y axes are the major and minor axes, respectively, and z is the orbital axis ± to 
the orbital plan e. Leng t hs are in semi-major axes and grayscale indicates surface 
brightness. See iMadural i20 1 Ol) for modeling details. 



To investigate the bore-hole effect, we use 3-D SPH simulations of t] Car's col- 
liding winds, together with a modified version of the visualization program SPLASH 
dPrice 2007) . to generate renderings of the surface brightness of ija in each of the var- 
ious wavebands B through L as a function of orbital phase. Opacities and primary star 
source functions as a function of radius in each waveband were taken from the spheri- 
cally symmetric model of ija by HOI, H06. Details can be found in Madura (2010) . 

Our models reveal that when tjb is at or near apastron, the WWIR is indeed too far 
from T]A and there is no bore hole (Fig. |2j left panel). As tjb moves closer to tja during 
its orbit, the WWIR gradually penetrates into the primary photosphere, creating a bore- 
hole effect that increases up until periastron passage (Fig. |2j middle panel), at which 
point ijB quickly wraps around the back side of ija, the bore hole briefly vanishing as it 
faces away from the observer. After periastron, the bore hole reappears on the opposite 
side of ijA (Fig.|2j right panel) and then slowly fades as 775 moves back toward apastron. 

Using the above 3-D model, we also generate synthetic photometric light curves in 



each waveband for comparison to the available ground-based photometry (IWhitelock et al. 
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Figure 3. Synthetic light curves from the bore-hole model (circles) com- 
pared to the photometric observations of rj Car's 2003.5 event (squares, 
iFernandez-Laius et al.ll2010h in the B (left) and K (right) bands . 



20041 : iFemandez-Laius etaDl201(]|) . We find that our bore-hole model reproduces the 



steep rise and drop observed in each waveband before periastron, and gives roughly 
the same peak-to-peak change in magnitude and duration of the 'eclipse-like' events 
seen every 5.54-years (see Figure [3] and Madura 2010). We thus find that t]b and the 
WWIR have a profound influence on t/^'s extended photosphere and the photometric 
observables. The bore-hole effect likely explains t] Car's periodic eclipse-like events. 



3. Effects of a Bore Hole on i] Car's Spectrum and Interferometric Observables 



Using the 2-D code of iBusche & HiUied (l2005l) . we have also investigated how a bore- 



hole effect could alter the observed optical and UV spectra of rj Car. Details on the 
modeling approach can be found in Groh (20 11) and references therein. 

Assuming a single-star scenario, HOI obtained a reasonable fit to the observed HST 
optical spectrum obtained just after periastron, (p ~ 0.05. The model spectrum from 
HOI reproduces well the emission line profiles of H, Fe II, and N II lines. However, 
compared to the observations, the HOI model overestimates the amount of P-Cygni 
absorption. The comparison is even worse as one moves toward apastron, when the 
observations show little or no P-Cygni absorption in H and Fe II lines. 

Using our 2-D radiative transfer model, which takes into account the cavity in 
the wind of tja caused by t]b, we computed the synthetic optical spectrum of rj Car at 
<p ~ 0.6 (see figure 1 of Groh 2011). We assume the same parameters for tja as HOI and 
a standard geometry of the WWIR at apastron as predicted by the 3-D SPH simulation. 
For a viewing angle with inclination / = 41° and longitude of periastron a> = 270°, the 
2-D model produces a much better fit to the P-Cygni absorption line profiles of H and 
Fe II lines than the 1-D HOI model, while still fitting the emission line profiles. The 
improved fit to the P-Cygni absorption is due to the cavity in the wind of t]a, which 
reduces the H and Fe II optical depths in line-of-sight to the primary. We find similar 
results when modeling rj Car's observed UV spectrum. 

A preliminary study of synthetic Ha line profiles generated using our 2-D cav- 
ity model of t]a also reveals broader, deeper P-Cygni absorption near the stellar poles. 
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Figure 4. Synthetic line profiles of Ha generated using the 2-D model discussed 
in the text, emphasizing the P-Cygni absorption at three stellar latitudes for tia- 



compared to the equator (Figure |4ll. A bore-hole or wind-cavity model may therefore 
provide an alternative explanation for HST long-slit spectral observations of Ha line 
profiles as a function of stellar latitude that ha ve been historical ly interpreted as evi- 
dence that the wind of tja is latitude dep endent (Smith et al.ll2003h . 

Finally, using a similar 2-D model, iGroh et all ( 20101 ) find that the density struc- 
ture of T/zi's wind can be sufficiently disturbed by tjb and the associated wind cavity, thus 



mimicking the effects of fast rotation in the interferometric observables. Gro h et al. 



do 10) further show that even if tja is a rapid rotator, models of the interferometric data 
are not unique, and both prolate- and oblate-wind models can reproduce current inter- 
ferometric observations. These prolate- and oblate-wind models additionally suggest 
that the rotation axis of ija would not be aligned with the Homunculus polar axis. 



4. Conclusions 

Our conclusions are rather straightforward; based on the results of our 3-D hydrody- 
namical simulations and multi-D radiative transfer modelling, we find that the sec- 
ondary star T]B and WWIR significantly affect the extended wind photosphere of tja, 
and thus multi-wavelength observations of the rj Car system. 

Acknowledgments. We thank Nathan Smith for discussions that led to this work. 
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